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ABSTRACT: We calculate the partition function of partially quenched chiral perturbation theory 
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1. Introduction 

At low energies, QCD can be described by a chiral effective theory. If the theory is considered 
in a finite volume and for small quark masses, the standard p-regime power counting is replaced 
by the e-regime power counting introduced by Gasser and Leutwyler [1]. To leading order in the 
e-regime, the partition function is dominated by the contribution of the zero-momentum modes 
of the Nambu-Goldstone (NG) bosons [1,2]. In this limit the theory becomes zero-dimensional 
and is therefore described by chiral random matrix theory (RMT) [3], see [4, 5] for reviews. The 
low-energy constants (LEC) appearing in the chiral effective Lagrangian, which are of great phe- 
nomenological importance, can be determined by fitting analytical results from RMT to lattice data 
for the eigenvalue spectrum of the Dirac operator. The lowest-order LECs are S and F. While 
X can be determined rather easily, e.g., from the distribution of the small Dirac eigenvalues, the 
extraction of F is somewhat more complicated and requires the inclusion of a suitable chemical 
potential [6, 7]. 

Since lattice simulations are restricted to a finite volume, it is important to take into account 
finite-volume corrections to the RMT results, which can be obtained by going to next-to-leading 
order (NLO) in the e-regime. Recently, finite-volume corrections to the unquenched partition func- 
tion of QCD in the e-regime were obtained in [8, 9]. However, in order to extract the relevant 
eigenvalue correlation functions the partially quenched partition function of QCD is needed. A 
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relatively simple method to obtain the partially quenched theory is to introduce n replicated flavors 
in the unquenched theory and then to analytically continue in the discrete number of quark flavors 
to zero. This so-called replica trick was first used in the theory of disordered systems [10]. It is 
potentially problematic since the analytic continuation from an isolated set of points is not uniquely 
defined. Nevertheless, a number of authors have succeeded to construct proper analytic continua- 
tions from which correct results could be obtained, see, e.g., [11, 12, 13]. Several publications in 
QCD have used the replica trick for perturbative calculations while borrowing exact result for the 
non-perturbative part of the theory from RMT [14, 15, 16, 17]. 

In this publication we choose to use an alternative way to obtain the partially quenched theory 
that does not suffer from the potential problems of the replica trick and can therefore be used 
to check and extend previous results. In addition to the sea quarks, we introduce fermionic and 
bosonic valence quarks. In nuclear physics and condensed matter physics this method is known as 
the supersymmetry method or Efetov method for quenched disorder [18]. In the context of QCD 
this idea was first used by Morel [19]. The effective low-energy theory of QCD with Nf + N v 
quarks and N v bosonic quarks was developed by Bernard and Golterman [20] and by Sharpe and 
Shoresh [21]. In this work we use the effective theory without a singlet particle as discussed by 
Sharpe and Shoresh and consider it in a finite volume and for small quark masses. In order to access 
F in addition to E, we include an imaginary quark chemical potential \i [6, 7]. (A first exploratory 
lattice study of this idea was performed in Ref. [22].) We compute the partition function at next- 
to-leading order in the e-regime and thereby obtain finite-volume corrections of order 1 / \/V to the 
partially quenched theory that translate into finite-volume corrections to the LECs E and F. Our 
results agree with previous results for the unquenched partition function [1, 8, 9]. As a side result 
we demonstrate that the parametrization of the NG manifold by Sharpe and Shoresh leads to the 
correct universal limit, in analogy to the results of Refs. [23, 24] where a different parametrization 
was used. 

An important question is to what extent the finite-volume effects in the determination of a par- 
ticular quantity, such as E or F, are universal in the sense that different methods used to determine 
this quantity give rise to the same finite-volume effects. In general the effects of the finite volume 
depend on the method, see, e.g., the finite-volume effects in the determination of F in Ref. [25]. In 
the present paper we show that at next-to-leading order in the e-expansion the partially quenched 
partition function is equal to its infinite-volume counterpart with E and F replaced by effective 
values E e ff and F e ff. Since the knowledge of the analytic form of the partially quenched partition 
function suffices to determine all spectral correlation functions of the Dirac operator Ip we find that 
all quantities that can be expressed in terms of spectral correlation functions of ij) give rise to the 
same finite-volume corrections to E and F. 

This paper is structured as follows. In section 2 we review the partially quenched theory and 
how it can be used to compute spectral correlation functions. We also review the corresponding 
effective low-energy theory in the formulation of Sharpe and Shoresh, both at fixed vacuum angle 
9 and at fixed topology v. In section 3 we compute the finite- volume corrections of order 1 / \/V to 
the partially quenched theory, and thus to E and F. We also show that the correct universal limit 
is obtained from the formulation of Sharpe and Shoresh. Conclusions are drawn in section 4. An 
appendix is provided to collect some useful formulas for the massless propagator in dimensional 
regularization, including commonly used shape coefficients. 
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2. QCD with Nf + N v quarks and N v bosonic quarks in a finite volume 



In this section we consider QCD with Nf + N v quarks and N v bosonic quarks (Morel's bosonic 
spin-1/2 ghost fields [19]) in a box of volume V = m ^ e Euclidean formalism. The 

temporal extent of the box is given by Lq, and thus the temperature of the system is T = 1/Lq. 
Unless stated otherwise we consider the partially quenched case of Nf > 0. 

2.1 The partition function and spectral correlation functions 

We define QCD with Nf + N v quarks and N v bosonic quarks by the partition function 



Z 



J d[A] 



-Sym 



■N f 

Jdet(0 + m/) 

./=i 



yr detjlp + m vi ) 
l = ldet(^ + m'J 



(2.1) 



where the integral is over all gauge fields A, Sym is the Yang-Mills action, Ip is the Dirac operator, 
mi, . . . , mjv- are the masses of the sea quarks, m v \, . . . , m V N v are the masses of the fermionic 



valence quarks, and m' vl , 



, m 



vN, 



are the masses of the bosonic valence quarks. By setting the 



mass m v i of a valence quark equal to the mass m' vi of the corresponding bosonic quark, the ratio 
of determinants of this pair cancels and the flavor i is quenched. 

Next we rewrite the determinants in terms of fermionic quark fields ip and bosonic quark fields 
ip using 



det(# + 



in) 



1 



and 



det(# + 



m) 



d[iptp] e' 



- f d 4 xip(Ip+m)ip 



J d 4 x ip(Ip+m)ip 



(2.2) 



(2.3) 



where tp and -ip are independent Grassmann variables with Berezin integral J d[V>V]> an< i <P an d 
(p are commuting complex fields related by complex conjugation, (p = ip*. The integrals in the 
exponents are over space-time. Note that the right-hand side of Eq. (2.3) only converges if all 
eigenvalues of Ip + m have a positive real part. Since Ip is anti-Hermitian this condition is satisfied 
as long as Re m > 0. Thus 



Z= d[A] d[*V] e -SYu-fd^9mM)-* 



(2.4) 



with mass matrix M = diag(mi , . . . , mjv, , m v i , . . . , m v N v , 



vli • 



,m: 



vN v 



) and fields 



ip <p 



(2.5) 



At nonzero temperature we have to choose anti-periodic boundary conditions in the temporal direc- 
tion for the fermionic quarks. A pair of fermionic and bosonic quarks at equal mass has to cancel 
in Eq. (2.1), and therefore we have to choose anti-periodic boundary conditions in the temporal 
direction also for the bosonic quarks (in the same way as Faddeev-Popov ghosts acquire periodic 
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boundary conditions at nonzero temperature [26]). This will amount to periodic boundary condi- 
tions for pseudo-NG fermions composed of quarks and bosonic anti-quarks (or of anti-quarks and 
bosonic quarks). 

The vacuum expectation value of an operator O is given by 

(0) = I J d [A] d[<2V] O e~ s ™~f d ' x . (2.6) 

For example, choosing N v = 1, the presence of a bosonic quark can be used to obtain the spectral 
density (or one-point function) of the Dirac operator 0, 

p(A) = (Tr S(lp-iX)) = lim-Re(Tr( J 0-iA + e)^ 1 ) , (2.7) 

£— >0 7T 



by using 



d 

(Tr(0 + m) _1 ) = — -logZ(mi, . . . ,m Nf ;m v ,m' v ) 



(2.8) 



m v =m'=m 



Analogously, higher-order spectral correlation functions can be obtained using N v = k, where 
k is the desired order. From these /c -point functions we can also compute individual eigenvalue 
distributions [27]. 

2.2 The effective low-energy theory at fixed vacuum angle 

In this section we briefly discuss how to determine the relevant low-energy degrees of freedom for 
QCD with Nf + N v quarks and N v bosonic quarks. For details we refer to Ref. [21]. The general 
procedure is as follows. We first determine the non-anomalous symmetries of the Lagrangian that 
act non-trivially on the vacuum. Then we restrict the remaining symmetry generators to a subset 
that is sufficient to generate all Ward identities associated with the flavor symmetries. This subset of 
symmetry generators then determines the relevant NG manifold of the effective low-energy theory. 
The Lagrangian of the quark sector is given by 

£ Q = *($ + M)*, (2.9) 

which in the massless case (M = 0) has vector and axial symmetries. The vacuum of this theory 
is invariant under the vector symmetry. The axial symmetry, however, acts non-trivially on the 
vacuum. 1 The axial symmetry is defined by a supermanifold [28] with base 

Gl(iV> + ^ 1 ,)®[Gl(JV t ,)/U(JV 1 ,)], (2.10) 

where Gl is the general linear group and U is its unitary subgroup. The factor Gl(iV/ + N v ) 
acts on the quark sector while Gl(N v ) / XJ (N v ) acts on the bosonic quark sector [23, 24]. The 
reason for the smaller symmetry group of the bosonic quark sector is that ip and Cp are related by 
complex conjugation, while ip and ip are independent in the functional integral. The measure of the 



1 For the detailed arguments concerning the symmetry breaking pattern of QCD with Nf + N v quarks and N v bosonic 
quarks we again refer to Ref. [21]. 
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functional integral restricted to the topological sector v transforms under axial transformations Ua 
as [29] 

-> Sdet iy (C/ A ) , (2.11) 

where Sdet is the superdeterminant [18]. Thus, for v ^ 0, only axial transformations with 
Sdet(LOi) = 1 leave the measure invariant, i.e., are non-anomalous. Let us express an arbitrary 
axial transformation U A by 

U A = exp (iG A ) = expi \ UA K ] , (2.12) 

where k and R are independent N v x (Nf + N v ) matrices with elements in the Grassmann algebra, 
ua lives in the group algebra of Gl(Nf+N v ), and u' A lives in the group algebra of Gl(A^) / \J(N V ). 
The restriction Sdet([/yi) = 1 amounts to the requirement of a vanishing supertrace [18] of Ga, 
i.e., Str Ga = Tr ua — TV u' A = 0. Next we restrict the remaining axial symmetries to the minimal 
subset that is necessary to generate all Ward identities of the full symmetry. Note that Gl(iVj + N v ) 
contains the same generators as U(Nf + 2V„) with real coordinates replaced by complex ones. 
Since this does not give rise to additional Ward identities it is sufficient to keep either the real or 
the imaginary part of each coordinate. The choice made in Ref. [21] is 

G = (* RT \ + iLf ( N '° 1n J+ n - | (2 13) 

A \Km'J y/(N f + N v )N v Nf \ (N f + N v )l Nv )' 

where ir = tt* and ir' = tt'' are traceless Hermitian matrices of dimension Nf + N v and N v , 
respectively, ip G R, and l n is the ?i-dimensional identity matrix. This choice leads to the correct 
signs of the kinetic terms of the NG particles in the effective low-energy theory and will also be 
used in the rest of this paper. Note that for Nf = also the flavor singlet particle will give rise to 
long-range correlations [21] and thus has to be included in the effective theory. 

The transformation properties of the massive theory under axial transformations as well as the 
Lorentz group now dictate the form of the Lagrangian of the effective theory [29] . To leading order 
in U(x), dpU(x), and M we find 

F 2 S 
£ eff = _ Str [dpUix^dpUix)} - - Str [M*U(x) + U(x)~ l M] , (2.14) 

where F and E are low-energy constants and the NG manifold U(x) is obtained by promoting the 
coordinates ir, tt', k, R, and ip in Eq. (2.13) to fields with 

U(x) = exp(iG A (x)) . (2.15) 

The theory in a 6- vacuum is then obtained by rotating the sea quark masses, 

£ eff (0) = £_ Str [dpU{x)- l d p U{x)\ - | Str [M^ e~ iS/N f U{x) + U {x)' 1 e iB 1 N f M] , (2.16) 



where 

" fl »> °\ (2.17) 
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is an (Nf + 2N V ) -dimensional matrix that projects onto the sea-quark sector. The partition function 
of the effective theory at fixed 9 is thus given by 

Z eS (6) = J d[U] e -/^£effW , (2.18) 

where d[U] is the invariant integration measure associated with the supermanifold [28]. We restrict 
ourselves to the effective theory in the rest of this paper and thus drop the subscript in the following. 

2.3 The effective low-energy theory at fixed topology u 

The partition function at fixed #-angle is given by the Fourier series 

oo 

Z(8)= e ie "Z„, (2.19) 

v= — oo 

and thus the partition function at fixed topological charge v is obtained by the Fourier transform 

Z v = — [ d9 e~ i6v Z[9) . (2.20) 
2tt Jo 

For the partition function defined in Eq. (2.18) this means 

Z v = J d9 J d[U] exp i-iOv - J d A x Str [d P U {x)~ l d P U {x)] 

- | Str [M^e~ iS ' N fU{x) + U(x)- x e a ' N 'M]\\ . (2.21) 

If we separate the constant mode Uq from U{x) by the ansatz 

U(x) = U exp(iG A {x)) (2.22) 
with J d 4 x Ga{x) = and Uq = exp(iG^), we can absorb 9 in Uo by 

9 ( lN f °\ n ™ 

vr ^ vro = vro - — Q ' Q j , (2.23) 

where itq is the constant mode of the pion fields in the fermionic quark sector of G\. To avoid 
confusion with (2.17) we mention that the matrix in (2.23) has dimension Nf + N v . Note that we 
absorb the #-angle only in the sea sector of the theory. This yields 

Z v = J d[U] Sdet^C/o) exp j- j d 4 x Str [d p U(x)- l d p U{x)] 

- | Str [M^U{x) + U{x)- 1 M]^j \ , (2.24) 

where the integration manifold for the constant mode is changed from (2.13) to 

G o = (*o R l\ , ifo ( N v ~b-N f +N v \ 25 

A UoiK) VWf + N v )N v Nf I (N f + N v )l Nv r 
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in which fro now generates \J(Nf + N v ) instead of SU(A r j + N v ) 2 while tt' , Rq, kq and ipo are 
defined in the same way as their counterparts in Eq. (2.13). Note that this parametrization of the 
constant mode is different from the parametrization used previously in the literature [23, 24]. In 
section 3.5 we will show that this parametrization again yields the universal RMT result. 

3. Finite-volume corrections 

3.1 The e-expansion in the effective theory with imaginary chemical potential 

For convenience we redefine the NG manifold with a different normalization of the fields by 

U(x) = U exp^j^€{x)J (3.1) 

with 

£( x \ = ( 7r ( x ) rT ( x )\ , Mx) ( N v t Nf +N v \ 

\k(x) iTr>(x)) y/(N f + N v )N v Nj { (N f + N v )t Nv J- 

The constant mode is separated in Uq, and thus f d 4 x^(x) = 0. For nonzero imaginary chemical 
potential the Lagrangian of the effective theory is given by 

F 2 £ 
C = — Str \y p U{x)- x V p U{x)] - - Str [M*U{x) + U(x)~ l M] (3.3) 

with 

V p U{x) = d p U(x) - iS p0 [C, U(x)} , (3.4) 

where C = diag(/^i, . . . , [j,N f , fj, vl , . . . , HvN v ,tf v i, ■■■ , /4iV"J and *M* is tne imaginary chemical 
potential of quark flavor i. We use the e-regime power counting [1] defined by 

V ~ e" 4 , M ~ e 4 , /i - e 2 , <9 p ~ £ , £(z) ~ e . (3.5) 

Note that the expansion in e 2 amounts to an expansion in 1/yV. To leading order in e 2 the 
Lagrangian is given by 

Co = \ Str \d p i{x)d p i{x)\ - | Str [M*U + U^M] - ^ Str [C, U^ 1 } [C, U ] . (3.6) 
The next-to-leading order terms in e 2 are 

C 2 = C% + C% + C% (3.7) 



2 The addition of 1 n } to the generators of SU(7V/ + N v ) suffices to generate U(JV/ + N v ). The normalization of 9 
in Eq. (2.23) yields the correct integration domain. 
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with 

£¥ = Str i Miu o^) 2 + t(x) 2 U l M] , (3.8) 

£ = -IstrU^CUol^x),^,^))) -^Str^CC/o + C)^),^)] , (3.9) 

$ = 7^2 Str[a^(x),^)][^(x),e(x)] - ^^Str ^[C, C/ ] [£(*), [d £(x), £(*)]] . 

(3.10) 

In this section we will integrate out the fluctuations in £ in order to obtain an effective finite-volume 
partition function. The term C% couples to Uq and M, and thus corrects the leading-order mass 
term. In section 3.3 we discuss its effect on the low-energy constant E. The term couples to Uq 
and C and corrects the leading-order chemical potential term. Its effect on the low-energy constant 
F is discussed in section 3.4. The first term in £^ can be ignored since it does not couple to Uq and 
therefore only amounts to an overall factor in the effective finite-volume partition function. The 
second term in £^ can be ignored at the order at which we are working since it does not give rise 
to leading-order corrections to £ or F. 

The integration measure for the parametrization of Eq. (3.1) is of the form 

d[U] = d[U }d[Z]J(0, (3-11) 

where d[Uo] is the invariant measure for the constant-mode integral, d[£] is the flat path integral 
measure of the fields £, and J{£) is the Jacobian corresponding to the change of variables of 
Eq. (3.1). Since £ does not contain constant modes the kinetic term in Eq. (3.6) suppresses large 
fluctuations in £, and thus the integrand vanishes at the integration boundaries of the it- and tt'- 
fields. Therefore the invariant integration measure is well-defined and there are no anomalous 
contributions by Efetov-Wegner terms [30, 31]. The Jacobian must be of the form 

J{i) = l + 0{e 2 ) (3.12) 

since there can be no contribution from a linear term in £ because of f d 4 x £(x) = 0. Thus, at 
next-to-leading order the Jacobian only contributes an overall factor to the effective finite-volume 
partition function. 3 

3.2 The propagator 

The kinetic term of the Lagrangian in terms of the fields it, it' , ip, R, and k is given by 

i Str [(d p 0(d p 0] = \ Tr [(d p 7r)(d p n)] + \ Tr [(V)(^')] + \ ^ [(d„<p)(d p <p)] 

+ {dpR j i){dpK ji ) . (3.13) 

Since the mass term C^ 1 °f the Lagrangian, see (3.8), is of order 0(e 2 ), the fields are effectively 
massless. The massless propagator without zero modes, A(x), is finite in dimensional regulariza- 
tion [32]. In appendix A we give explicit expressions for the relevant propagators used in this work. 



3 At higher orders in e the effects of the Jacobian can no longer be absorbed in an overall prefactor of the partition 
function. 
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For the pion fields tt and n' the propagators are given by [33, 9] 



(vr(x) afe 7r(y) cd ) = A{x - y) 
(n'{x) ab 7r'(y) cd ) = A{x - y) 
where the average is defined by 



bad&bc 



1 



N f + N v 
1 



babied 



N, 



SnhS, 



aWcd 



fd[£\ Oft] e~f d4xC ° 



For the scalar field tp and for the fermionic field k the propagators are easily shown to be 

(R(x) ab K(y) cd ) = -A(x - y)5 ac S bd , 
(<p(x)<p(y)) = A(x - y) . 



Using the identities 



i 



+ 



Nf, 



1 



N f + N v {N f + N v )N f N v N f ' 

l_ {Nf + N v ) 2 J_ 

N v + (N f + N v )NfN v N f ' 

we thus find the propagator of the composite field £ to be 



(£{x)ab€(y)cd) = A(x - y) 



1 



S a d8bc(-i) Eb - -ir^abScd 
Nf 



with 



(3.14) 
(3.15) 



(3.16) 



(3.17) 
(3.18) 



(3.19) 
(3.20) 



(3.21) 



(3.22) 



for 1 < b < N f + N v , 

1 for Nf + N v <b<Nf + 2N V . 

Note that there is no explicit dependence on the number N v of valence quarks in this propagator. 
3.3 Finite-volume corrections to £ 

We now integrate out the fluctuations in the 0{e 2 ) mass term C^ 1 to obtain the finite-volume 
corrections to the leading-order mass term in Cq. Using (3.21) it is straightforward to show that 



(Str[AZ(x)BZ(y)]) = A(x-y) 

By expanding the action we find that the term 

E 



StrAStrS 



1 



Str AB 



(3.23) 



J d^x L^^[M^(xf + ^xfU^M) 



(3.24) 
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corrects the leading-order mass term in the Lagrangian, 

- - Str [M*U + U^M] , (3.25) 

to 



£ 
2 



N 2 -l 



Str [M*U + Uq X M] . (3.26) 



Thus at next-to-leading order we can read off an effective low-energy constant S e ff given by 

T « Nr - 1 _ 

-W = 1 - A^ A(0) • (3 " 27) 

This is the same result as previously derived for the unquenched partition function [1, 8, 9]. It can 
be shown that at next-to-next-to-leading order (NNLO) it is no longer possible to absorb the effects 
of the finite volume in an effective low-energy constant £ e ff. 

3.4 Finite-volume corrections to F 

The calculation of the finite-volume corrections to F is slightly more involved. The non-vanishing 
corrections to the leading-order imaginary chemical potential term are given by Eq. (3.9). We first 
calculate the contribution of the first term in (3.9), 

-\ j d 4 x (Str U^CUoUx),[C,^x)}}) 

= -\jd 4 x (Str U Q l CU [2^x)C^x) - £(x) 2 C - C£(x) 2 ]) 

= -yA(0) [(Str C) 2 - N f Str U^CUqC] , (3.28) 

where we have used (3.23). The first term in (3.28) couples only to C 2 and thus amounts only 
to a prefactor in the effective finite-volume partition function. The correction to the leading-order 
Lagrangian obtained from (3.28) is thus given by 

^1 N f Str [C,U^][C,U ]. (3.29) 
The contribution of the second term in (3.9) is given by 



J d 4 x (Str (U^CUo + C)[t{x), d £(x)}) ~ S A(0) = (3.30) 



due to the symmetry A(x) = A(— x). However, the square of this term gives a nonzero contribu- 
tion. We need to calculate 

d 4 x Str([/ - 1 CC/ + C)[e(x),^(x)] V 
1 



2 \ V 2 



d 4 x / A(Str(y[e(x),5 e(x)])Str(yK(y),5 e(y)])> (3.31) 
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with Y = U 1 CUq + C. After performing all relevant contractions using (3.21) we find 



(Str[Yt(x)t(x')] Str[^(y)C(2/0 



'o 



1 



A{x - x')A{y - y') |@tr YfNj - 2(Str Y) 2 + ^(Str Y) 2 
+ A(x-y)A(x'-y') [(StrY) 2 



N 



StrY 2 + ^(StrY)' 



A(x - y')A(x - y) 



N f Str Y 2 



N 



StrY 2 + 



N2 f 



(StrY) 2 



(3.32) 



Since Str Y = 2 Str C does not couple to Uq we only need to take into account the terms involving 
Str Y 2 . We denote the irrelevant terms by ". . ." and write 



(3.32) 



StrY 2 



N f 

We need to calculate 



[2A(x - y)A(x' - y') + (2 - N 2 f )A(x - y')A(x' -y)] + ... (3.33) 



(d x , Q - d X0 )(d y , - d yo ) (Stv[Y^x)^')] Str[Ye(y)e(y')]> c 



x=x' , y=y' 



-2N f StrY 2 [(d A(x - y)) {d A(x - y)) - (a o 2 A(x - y))A{x - y)] + . . . (3.34) 



Thus we find 



(3.31) = ~N f StrY 2 J d 4 x (d A(x)) 2 + 



(3.35) 



where we have used the fact that the propagator is periodic in time. Therefore the corrections to 
the effective Lagrangian are given by 



^NfStriCU^iCUo] J d 4 x {d A(x)Y 



(3.36) 



Combining (3.29) and (3.36), we find that the fluctuations correct the leading-order contribution to 
the Lagrangian, 



StrlCU^mUo] 



(3.37) 



to 



Stv [C,U^}[C,U ] 



2Nf 
F 2 



A(0) - J d 4 x(d A(x))'- 



Thus at next-to-leading order we find an effective low-energy constant F e g given by 



-Feff 



N f 



1-^(A(0)- I d 4 x{d A(x)) 



(3.38) 



(3.39) 



This again agrees with the result for the unquenched partition function [8, 9]. As in the case of S, at 
NNLO it is no longer possible to absorb the effects of the finite volume in an effective low-energy 
constant F e ff. 
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3.5 The universal limit 

In this section we concern ourselves with the limit V — > oo while keeping MVH ~ O(e ). It is 
well known that QCD in this limit behaves in a universal way and agrees with chiral RMT. It was 
first shown in Refs. [23, 24] how universal results for the Dirac spectrum (obtained earlier in RMT) 
can be derived from the effective low-energy theory. In the following we show that the correct 
universal limit also follows from the effective theory in the formulation of Sharpe and Shoresh 
described above. For simplicity we restrict ourselves to the case of vanishing imaginary chemical 
potential, C = 0. 

Since the fluctuations in £ are suppressed for V — > oo only the zero-mode integral survives in 
this limit, and the partition function for fixed topological charge u is given by 

Z u = J d[U ] Sdef{U ) exp Str [M*U + Uq X M] \ , (3.40) 

where the integration manifold is specified in Eq. (2.25). There are different methods to calculate 
integrals over supermanifolds, see, e.g., [34, 35, 36]. In our case it is sufficient to choose an 
explicit parametrization and reduce the integral to ordinary group integrals. For convenience we 
use a slightly different notation and calculate 

Z„ = j d[U] Sdet"(t0 exp (str [M*U + ET^M]) (3.41) 

with integration manifold given by 

(Ve Nvt f \ fo R T \ 
U ={ V.W«M») a *(*o)= lW '" <3 " 42) 

where V G U(JV> + N v ), V G Gl(N v ) / U(N V ) with det V = 1, and ip G R. Thus we have 

Sdet U = det V = e id (3.43) 

with 9 G [0, 2tt). This is the zero-mode integral following from the parametrization used in the 
perturbative calculation above. In the literature a similar integral was computed to determine the 
static limit of partially quenched chiral perturbation theory [23, 24] that amounts to replacing U c 
by 

u < - (o v**) ' < 3 44) 

Note first that a parametrization such as U = U c U g above leads to factorization of the correspond- 
ing measure as 

d[U] = d[U c ]d[U g ] . (3.45) 

This is due to the fact that the invariant length element is 

ds 2 = Str [dUdiU' 1 )} 

= Str [dUcdiU' 1 ) + dUgdiU' 1 ) - 2U~ 1 dU c dU g U- 1 } 

= Str [dU c d{U~ l ) + dUgdiUg 1 )} (3.46) 
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since 



- (I f ) • <" 7 > 

U~ l dU c is block diagonal, and therefore 

StrfC/- 1 dUcdUgUg 1 } = . (3.48) 
In both parametrizations the measure of V, V', and <p also factorizes. Thus 

d[U] = d[U g ]d[V]d[V'}d(p (3.49) 

in both cases. Note that this parametrization has no contributions from Efetov-Wegner terms, as 
was discussed in a special case in the literature [24]. Introducing the short 4iand notation 



we find for the first parametrization 

Str [M ] U + MU^ 1 } = Str [M ] U c U g + MU^U' 1 } = Str [U C X + U^Y] 

= Tr [Ve N ^X s - V'e {N f +N ^X bb + V^e'^Yn - V^ +J ^*Ybb] . (3.51) 

Next we use a result of [37], 

[ d[V] det"(C0 exp \Tv(AU + BU' 1 )) = c p det^" 1 )^ det N (3 . 52) 
7u(p) A (M ) 

where c p is a constant, A(// 2 ) is the Vandermonde determinant, and the /x| are the eigenvalues of 
AB. Thus the integral over V results in 



-N v {N f +N v )u V(]pt(YttX -i Y /2 det[^ 1 J y+J -_i(2/^)] 



e-^^+^^detCYffX- 1 )^ 2 1 ? A 7,T 2 r (3 - 53) 
with /i 2 the eigenvalues of X^Y^. In the second parametrization we find 

Str [M*U + Ml/- 1 } = Tr [VX S - V'e v/Nv X hh + F" 1 ^ - V'^e'^Y^} . (3.54) 

Note that in this parametrization we also have an additional factor of e~ vv from the superdetermi- 
nant. Thus the integral over V leads to 



-1^/2 det[^ 1 J y+J -_i(2/i i )] 
A(M 2 



e"^ detC^X- 1 )^ 2 J (3-55) 



with /i 2 already defined above. Now we let <p — > ipN v (Nf+N v ) in order to have the same prefactor 
of V' and F'" 1 in the supertrace. In both parametrizations the resulting integral is 

/ dp g ]d[V] r d^e-^ N f +N ^det{Y n X^y/ 2 ^r l(2/ii)) 

x exp (- TV [v'e^ N f +N ^X bb + V^"^^] ) . (3.56) 
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This completes the matching with Refs. [23, 24] and is sufficient to show that the parametrization 
of the NG manifold used in this work leads to the correct universal limit. 

In order to extend this proof to the general case of C / we would need to calculate the group 
integral 

/ d[U] det v (U) exp [Tr(AU + BU^ 1 ) + ^(DUDU^ 1 )] , (3.57) 

■/U(p) 

where A, B, and D are arbitrary complex p x p-matrices. This, however, is beyond the scope of 
this work. 



4. Conclusions 

In this work we have calculated the partially quenched partition function of QCD at next-to-leading 
order in the e-expansion at nonzero imaginary chemical potential. We considered a theory with 
Nf + N v fermionic quarks and N v bosonic quarks, as formulated by Sharpe and Shoresh [21], in 
a finite volume V with microscopic quark masses M, i.e., MVT, = O(e ). The knowledge of the 
analytic form of the partially quenched partition function suffices to obtain all spectral correlation 
functions of the Dirac operator J]). In this sense our results for the finite- volume behavior of the 
theory hold universally for all observables that can be obtained from spectral correlation functions 
of 0. We found that the partially quenched partition function has the same finite-volume correc- 
tions as the unquenched partition function of QCD with Nf quarks, i.e., at next-to-leading order in 
e there are effective low-energy constants £ e ff and F e ft, 

£ ft N 2 f - I _ 

-f = 1 " i^T*(0) , (3.27) 



_ n N f 

F 



1 - ^ (A(0) - J d 4 x {d A(x)) 2 ) , (3.39) 



where A (a;) is the massless propagator. In appendix A we give closed formulas for the relevant 
propagators in dimensional regularization and numerical values for typical geometries. As a side 
result of our calculation we showed that the constant-mode integral of this theory agrees with previ- 
ous results from random matrix theory. Therefore the correct universal limit is obtained at V — > oo. 
Note that the proof was only given for vanishing chemical potential and that the knowledge of the 
group integral (3.57) is needed to complete the proof also for nonzero chemical potential. 

In figure 1 we show the finite- volume corrections at NLO to the low-energy constants £ and F 
as a function of the box size L in a symmetric box. Note that the effects of the finite volume increase 
with the number of sea quark flavors Nf and that, depending on Nf, a box size of 3 — 5 fm is 
necessary to reduce the effects of the finite volume at NLO to about 10%. The effects are calculated 
at F = 90 MeV. In figure 2 we show the effect of an asymmetric box with Nf = 2 and L = 2 
fm. An important message of this figure is that the magnitude of the finite-volume corrections 
can be significantly reduced by choosing one large spatial dimension instead of a large temporal 
dimension. The reason for this behavior is that the chemical potential only affects the temporal 
direction, see Eq. (3.4), and therefore breaks the permutation symmetry of the four dimensions. 
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Figure 1: Volume-dependence at NLO of the low-energy constants £ e ff (left) an d Fes (right) m a symmetric 
box with dimensions Lq = L\ = L2 = £3 = L at F — 90 MeV. 
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- F etf /F 

— 2 eff /E 



1 








Figure 2: Effect of an asymmetric box with parameters Nf = 2, L = 2 fm, and F = 90 MeV. We 
compare a large temporal dimension Lq with L\ = L2 = L3 — L (left) to a large spatial dimension L3 with 
L = Lx = L 2 = L (right). 

This manifests itself in the propagator 

J d 4 x (d A(x)) 2 (4.1) 

which, as shown in Eq. (A.2), contains a term proportional to L^/y/V, where Lq is the size of 
the temporal dimension. This term leads to an enhancement of the corrections in case of a large 
temporal dimension. Choosing instead one large spatial dimension, the finite-volume corrections 
are reduced, unless the asymmetry is too large. For the parameters used in figure 2, the optimal 
value is L%/L w 2. 

This is good news. Many lattice simulations (at zero chemical potential) are performed with 
Li = L2 = L3 = L and Lq = 2L. To determine F, it suffices to introduce the imaginary 
chemical potential in the valence sector. Therefore, one can take a suitable set of existing dynamical 
configurations and redefine Lq <-> L3 before adding the chemical potential. 4 This will minimize 
the finite- volume corrections for both E and F, at least for the parameter values chosen in figure 2. 

4 Note that this procedure increases the temperature of the system by a factor of two. One needs to check that the 
system does not end up in the chirally restored phase, in which our results no longer apply. 
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A. The massless propagator in dimensional regularization 



For convenience we collect in this appendix explicit formulas for the massless propagator in di- 
mensional regularization, A(x), that were derived in Refs. [32, 33]. The two relevant quantities for 
the finite-volume corrections to E and F are given by 



A(0) 



ft 
V 



(A.1) 



and 



J d 4 x (d A(x)Y 



2W 



ft ~ ^kh 



V 



(A.2) 



where ft and fcoo are so-called shape coefficients, i.e., they only depend on the quantities Zj 
Li/V l / A with i = 0, 1, 2, 3. The shape coefficient ft is given by 



ft 



l r 

47T 



2 - - d-itfT 1 ) 



with 



dt f 



a-i(xj) 

where S(x) is an elliptic theta-function defined by 

oo 

S(x) = £ e 

n=—oo 

The shape coefficient fcoo is given by 

1 -E 

19 /-^ 



l[S(x 2 3 /t)-l 

.3=0 



COO 



12 4sinh(Z gn/2) 2 ' 

where the sum is over all integers (m, ^2,^3) 7^ (0, 0, 0) and 

3 

3=1 



(A3) 



(A.4) 



(A.5) 



(A.6) 



(A.7) 



In tables 1 and 2 we give numerical values for common shapes. 
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Lo/L 


1 


2 


3 


4 


Pi 


0.1404610 


0.0836011 


-0.0419417 


-0.215097 


&00 


0.0702305 


0.0833122 


0.0833333 


0.0833333 



Table 1: Coefficients for an asymmetric box with L\ = L-2 = L3 = L and temporal dimension Lq. 



L 3 /L 


1 


2 


3 


4 


Pi 


0.1404610 


0.0836011 


-0.0419417 


-0.215097 


koo 


0.0702305 


-0.0322630 


-0.2984300 


-0.731240 



Table 2: Coefficients for an asymmetric box with Lq = L\ = L2 = L and spatial dimension L3. Note that 
/?! is symmetric under the exchange of the temporal with a spatial dimension. 
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